Spin moment over 10-300 K and derealization of magnetic electrons above the 

Verwey transition in magnetite 
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In order to probe the magnetic ground state, we have carried out temperature dependent magnetic 
Compton scattering experiments on an oriented single crystal of magnetite (FeaCU), together with 
the corresponding first-principles band theory computations to gain insight into the measurements. 
An accurate value of the magnetic moment ns associated with unpaired spins is obtained directly 
over the temperature range of 10-300K. /is is found to be non-integral and to display an anomalous 
behavior with the direction of the external magnetic field near the Verwey transition. These results 
reveal how the magnetic properties enter the Verwey energy scale via spin-orbit coupling and the 
geometrical frustration of the spinel structure, even though the Curie temperature of magnetite is 
in excess of 800 K. The anisotropy of the magnetic Compton profiles increases through the Verwey 
temperature T v and indicates that magnetic electrons in the ground state of magnetite become 
delocalized on Fe B-sites above T v . 



T3 

C 

O 
o 



> 

in ■ 

o : 

m . 
m ■ 
O ■ 

: 

o . 

■ 



I 

c 

o 
o 



X 



PACS numbers: 73.22.Dj, 75.75.+a, 75.10.-b 

Magnetite was known to the ancient Greeks and its 
"magical" properties have been a source of fascination 
for millenia. In recent years, magnetite has become an 
archetype of complex correlated materials which display 
competing orders and interesting charge and orbital or- 
dering effects in a solid state setting. This complexity 
is accentuated by the expectation that Fe in the inverse 
spinel crystal structure of magnetite will assume a mixed 
valence. Simple chemical arguments suggest that on the 
tetrahedral sites Fe occurs as ferric cations (Fe 3+ ), while 
the octahedral sites are populated equally by ferric and 
ferrous (Fe 2+ ) cations. Also, magnetite is important for 
technological applications because it magnetizes sponta- 
neously like a ferromagnet but possesses low electrical 
conductivity so that it can be used at high frequencies 
where eddy currents make the use of iron impossible. 

Magnetite has also drawn considerable attention be- 
cause it displays a first order phase transition-the Ver- 
wey transition-at T V =120K, below which the resistivity 
increases by a factor of ~100 0,0. The earliest interpre- 
tation was that this is an order-disorder transition where 
the Fe 2+ and Fe 3+ cations on the octahedral sites become 
ordered below T v . However, this interpretation and var- 
ious pictures of electronic ordering on different Fe sites 
still are quite hotly debated issues 0, 0, H, H, • Band 
computations based on the local density approximation 
(LDA) show that the bonding electrons are shared in 
states that cannot be assigned solely to iron or oxygen 
Hi B EH and predict magnetite to be a half-metal above 
T v where the minority spin electrons are conducting, but 



the majority spins are insulating. 

It is clear that the behavior of magnetic Fe 3d elec- 
trons in magnetite is the key to understanding the prop- 
erties of this material. With this motivation, we have 
carried out a magnetic inelastic x-ray scattering investi- 
gation of an oriented single crystal of magnetite in the 
deeply inelastic (Compton) regime. In this way we obtain 
a direct measurement of the magnetic moment associated 
with unpaired spins in magnetite over the wide temper- 
ature range of 10—300 K. Notably, magnetic Compton 
scattering (MCS) is a genuine bulk probe which can de- 
tect the spin magnetic moment (ns) with a sensitivity of 
~0.1% (1 /Zfl per 1000 electrons) [llj. Previous estimates 
of ns are less direct, being based on the application of 
sum rules for extracting orbital moments (no) from x-ray 
magnetic circular dichroism (XMCD) spectra [l2| • More- 
over, a recent XMCD study [H} of magnetite reports a 
nearly vanishing no and an integral value of ns, which 
is in sharp contrast to the earlier XMCD results 12j of 
a large value of no and a non-integral ns 0,0. It is 
important therefore to determine independently whether 
ns is integral or not and no is small or large in order to 
understand the role of spin-orbit coupling in magnetite. 

An MCS experiment is well known to couple directly 
to the wavefunction of mag netic electrons in the ground 
state of the system 11 , 16[ 17| • The present MCS mea- 



surements over the 10-300K range thus allow us to com- 
pare the characteristics of the magnetic states in mag- 
netite above and below the Verwey transition. In particu- 
lar, we show that the anisotropy of the MCS spectrum be- 
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comes significantly larger as the system goes through T v . 
In order to help interpret our measurements, we report 
extensive computations of the magnetic momentum den- 
sity within the first-principles band theory framework. 
These and other model computations enable us to de- 
duce that the anisotropy of the MCS provides a handle 
for ascertaining the localized vs delocalized nature of the 
magnetic electrons in the ground state of magnetite. 

In an MCS experiment one measures the magnetic 
Compton profile (MCP), J m ag(Pz), for momentum trans- 
fer along the scattering vector p z , which is defined by 

Jmag(Pz) = Jt(Pz) ~ Jl(Pz) , (1) 

where J-j- (Jj) is the majority(minority) spin Compton 
profile. J ma g can be expressed in terms of a double inte- 
gral of the spin density, p mag (p)- 

Jmag (,Pz ) / / Pmag 

(p)dp x dp y , (2) 

where p m a g (p) = Pt(p) — Pl(p) is the difference of the 
majority and minority spin densities. We emphasize that 
the area under J m a g (Pz) gives the spin moment ps- In 
contrast, the information about the nature of the ground 
state wavefunction is described by the shape of J m ag{Pz)- 

Experiments were performed on bcamline 11-B at the 
Advanced Photon Source [3] using an elliptical multi- 
pole wiggler to generate circularly polarized photons. A 
high quality single crystal of Fe304 was used and MCPs 
were measured along the [100] and [110] crystal direc- 
tions at temperatures of 10K, 100K, 140K, and 300K. 
The Verwey transition temperature of the sample was 
determined via resistivity measurements to be 120 ± 1 
K. The incident photon energy was 125 keV. The scatter- 
ing angle 9 was 170°. All measurements were made under 
an external magnetic field of 7T oriented along the direc- 
tion of the scattering vector. The momentum resolution 
is estimated to be 0.4 a.u. (full- width- at-half-maximum). 
Magnetic calibration measurements [l8| were performed 
using a crystal of Fe oriented along [110]. 

The electronic structure was computed for the high 
temperature (Fd3m) phase of F23O4 within an all- 
electron, fully charge and spin self-consistent LDA-based 
band theory framework [191 ] . All computations were car- 
ried out to a high degree of accuracy, e.g., the crystal 
potential was converged to better than 10 -5 Ry. Our 
band structure is in good overall accord with the few 
available studies [1, HJ Various calculations agree 
on a half-metallic ferromagnetic band structure with the 
Fermi level crossing the minority spin ti g bands originat- 
ing from the octahedral Fe sites. The computed band 
structure and wavefunctions were used to obtain spin- 
resolved momentum densities p-f(P) an d Pi(p)- These 
quantities were computed over a fine mesh of 855 x 10 6 
points within a sphere of radius ~ 13 a.u. in p space 
[201 ] in order to properly account for the fine structure 
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FIG. 1: (Color online) J mag (Pz) for p z along [100] (blue dots) 
and [110] (red dots) as a function of temperature. Differ- 
ent datasets are offset vertically with respect to one another. 
Symbol size is representative of error bars. Inset gives the 
magnetic moment of unpaired spins fper formula unit) ob- 
tained from the area under the MCP [llfl . 



in p| and pi, and provided the dataset from which var- 
ious projections of the spin-resolved momentum density 
and the MCPs along the high symmetry directions were 
computed. 

The experimental spectra are summarized in Fig. [1] 
The area under the MCPs yields the magnetic mo- 
ment ps associated with uncompensated electronic spins 
shown in the inset, ps has a value of 3.54 /^/formula 
unit at 10K for the magnetic field along [100]. The un- 
certainty in ps is estimated to be ±0.05 ps- By compar- 
ing with the bulk magnetic moment of 4.05 pb [21] , this 
yields a corresponding orbital component of 0.51 ± 0.05 
Pb- Interestingly the 140K dataset near the Verwey 
transition displays a striking anomaly in that ps is 3.83 
Pb when the 7 T field is along [100], but for the field 
along [110] it is 3.28 pb- This directional preference of 
ps weakens again at higher temperatures. 

It is interesting to note that ps is non-integral. A sim- 
ple ionic picture gives ps — ^Pb per formula unit 0]. 
In general, band theory can accommodate non-integral 
ps since bands can be occupied partially. But, if spin is 
assumed to be a good quantum number, it is straight- 
forwardly argued that /ij in a half- metallic system must 
also be strictly integral. However, a non-integral ps value 
can be obtained if the spin-polarized band calculation is 
supplemented by spin-orbit interactions. The spin-orbit 
coupling has been estimated to be of the order of 10 meV 
[23 ] . so that we would expect a substantial orbital mo- 
ment due to Fe 2+ . 
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Bulk magnetization measurements in magnetite 23 



24j reveal a significant magneto-crystalline anisotropy. 
This anisotropy is relatively small above T v with [111] be- 
ing the easy axis, while below T v , the anisotropy is larger 
with [001] being the easy axis. However, the maximum 
saturation field is about 2T, so that under the field of 
7T used in the present measurements, we would not nor- 
mally expect anisotropic effects. Therefore, changes in 
the spin moment with the direction of the magnetic field 
seen in the inset to Fig. [TJ which are especially strik- 
ing at 140K, are anomalous and do not fit within the 
standard picture. The origin of this anomaly is unclear, 
but it may be related to the effect of geometrical frus- 
tration on magnetic properties above T v p}. Notably, 
Ref. [SH considers a Hamiltonian for the pyrochlore B 
lattice with a spin-orbit interaction. In this model, sev- 
eral complex orbital ordered states with noncollinear or- 
bital moments are found as a possible ground state for 
the high-temperature phase. This degeneracy of the low- 
est energy states can be broken by an external magnetic 
field in a manner which depends on the field direction. 
Since the energy differences between these states are very 
small, a magnetic field of 7T can be sufficient to switch 
between the states 26}. Further MCS measurements of 
the spin moment in magnetite over a range of temper- 
atures around T v should prove worthwhile in this con- 
nection. Here a note should also be made of the role 
of magnons, which will reduce the moment below the 
canonical value of 4 [is at high temperatures [271 ] . 

Some insight into the shape of the MCPs can be ob- 
tained through the LDA-based first-principles band the- 
ory computations. This is illustrated in Fig. [21 which 
shows a good level of agreement between the shapes of 
the computed and measured [110] MCPs. The LDA cor- 
rectly reproduces the overall width of the measured MCP 
as well as the pronounced peak at 1.3 a.u. The shallow 
bumps in the experimental MCP at momenta between 
2 — 3 a.u. and 4 — 5 a.u. are similar to features in the 
theoretical MCP. This fine structure is the hallmark of in- 
terference terms that cannot be explained within a purely 
atomic picture. 

In order to extract wavefunction localization features, 
we now discuss the anisotropy of the MCP defined as the 
difference between the two directional MCPs [28( , 
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FIG. 2: Shape of the 140K [110] experimental MCP of 
Fig. [1] is compared with corresponding theoretical (resolu- 
tion broadened) MCP obtained within the conventional LDA- 
based band theory framework. 



the Verwey transition, shows that although the LDA re- 
produces the overall features in the measured anisotropy, 
the amplitude of the anisotropy given by the LDA is too 
large by about a factor of 4 (note that plotted LDA curve 
is scaled down). In order to gain insight into the shape 
and amplitude of AJ magi we have considered this quan- 
tity for a variety of clusters of Fe A- and B-site atoms 
and their nearest neighbor (NN) Fe atoms, using linear 
combinations of Slater type orbitals (STOs) where the 
cluster wavefunction is constructed by allowing STOs on 
NN Fe atoms to mix via a small admixture parameter / 
[30I ] . Our results show that the shape of the experimental 
AJ mag , including the location of the peak at 1.12 a.u., 
can only be reproduced for a B-site Fe cluster in which 
the NN Fe atoms are situated at a NN distance dss=5.61 
a.u. (2.97 A) along the [110] direction [31|. In this case, 
the momentum wavefunction obtained by Fourier trans- 
forming the cluster wavefunction contains an oscillating 
factor of [1 + 2/ cos{cIbb(Px +Pt/)/v / 2}] 7 which yields the 
correct shape, and as Fig. [3] shows, for / = 0.1 the com- 
puted curve (dashed) gives a good agreement in shape 
and amplitude with the 140K and 300K measurements 



where the brackets 
perature 



denote an average over tern- in the low momentum region. We note that AJ mag is 



Jmag) is normalized at each temperature 



to the integral of the corresponding J, 



110] 

mag , 



AJmag integrates to zero. 



so that 

The average ( Jma°J ) pro- 



vides a temperature-independent "background" for high- 



[110] 



as a function of 



lighting changes in the shape of Jmag 
temperature. [29( 

Fig. [31 which focuses on the 140K and 300K data above 



n> ag 

rather insensitive to the inclusion of covalency with O 
atoms, and furthermore, the lack of oscillations in the 
cluster model at p z > 2 a.u. is related to inaccuracy in 
our STOs near the nucleus and is of little interest in the 
present context of solid state effects. 

In order to delineate the temperature dependence of 
AJ mag the results for the 10K and 100K datasets be- 
low T v are presented in Fig. 2] The peak located at 
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FIG. 3: Anisotropy of the MCP defined as AJ„ 



mag ^ — Jma°J , where the areas under all profiles are nor- 
malized to one. Experimental results for the f40K and 300K 
datasets are compared with the anisotropy obtained within 
the conventional LDA-based band theory framework and the 
model cluster computation discussed in the text. Note, the 
LDA curve has been scaled down by a factor of 4. 
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FIG. 4: Anisotropy of the MCP defined as AJ mag = 
Jmag , where the areas under all profiles are nor- 
malized to one. Experimental results for the 10K and 100K 
datasets are compared with the simple atomic analysis of 
Ref. [13]. The LDA curve is reproduced from Fig. [3] for refer- 
ence. 



1.12 a.u. is seen now to have essentially disappeared and 
the anisotropy is strongly reduced compared to even the 
scaled LDA curve. The rather flat result of Fig. |3] is con- 
sistent with a value of the admixture parameter i_~ 
[321 ] and with the simple atomic analysis of Ref. [13] ■ A 
similar result for / was obtained from an analysis of a 
positron annihilation study on magnetite carried out be- 
low T v [13]. 

To be more quantitative we define the anisotropy am- 
plitude 



A = 



\AJ mag (p)\ dp, 



(4) 



where p max — 5 a.u. is a momentum cutoff [33|. As we 
have already discussed in connection with Figs. [3] and [3] 
above, the strength of the peak at 1.12 a.u. in AJ mag (p) 
is related to the extent to which magnetic electrons on 
B-sites develop phase coherence via the admixture pa- 
rameter / on NN Fe B-sites. It may then be argued that 
the quantity A defined via the integral of Eq.[4]is a mea- 
sure of the area under the peak in AJ mag {p), and thus 
of the number of magnetic electrons participating in this 
coherence process [34 1. 

Fig. [5] shows that the value of A increases systemati- 
cally in going from the 10K to the 300K dataset. Inter- 
estingly, this behavior of A correlates with that of the ad- 
mixture parameter / discussed above in connection with 
Figs. |31 and 2] in that, like A, the parameter / also in- 
creases with temperature. In other words, the Verwey 



transition is accompanied by a change in the character 
of the ground state wavefunction such that the wavefunc- 
tion becomes significantly more delocalized on Fe B-sites 
above T v . Notably, the present analysis suggests that 
the B-site is involved in the change in valence across the 
Verwey transition as in Ref. [35(, and not the A-site as 
argued in Ref. |36] recently. Positron annihilation ex- 
periments [37I |38| . which also give information on spin- 
resolved momentum density, have revealed covalency ef- 
fects between Fe and O atoms but have failed to detect 
significant changes at T v . This result however may reflect 
the tendency of the positron to sample mostly O-ions and 
interstitial spaces rather than the Fe-ions. 

In conclusion, we have presented magnetic Compton 
scattering measurements on a magnetite single crystal, 
together with corresponding predictions of the MCP 
within the conventional band picture. An accurate value 
of the unpaired spin moment us is thus obtained directly 
over a wide temperature range. For example, at 10K /ig 
has clearly a non integer value of 3.54±0.05 /Lts/formula 
unit for the magnetic field along [100] demonstrating a 
non vanishing spin-orbit coupling. The ground state of 
magnetite is shown to be remarkably sensitive to the di- 
rection of the external magnetic field around the Ver- 
wey transition and to display an anomalous magnetic 
moment, which may be a manifestation of a large degen- 
eracy and associated geometrical frustration in the spinel 
lattice. The amplitude of the anisotropy of the MCP is 
shown to increase with temperature. We argue on this 
basis that the ground state wavefunction in the system 
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FIG. 5: Amplitude A of the MCP anisotropy obtained from 
Eqs. 3 and 4 as a function of temperature. 



becomes delocalized on Fc B-sites above T v . 
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